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FARADAY  ROTATION  NEAR  TEE  TRANSVERSE  REGION  OF  THE  lONOSroSRE 

ABSTRACT 

It  Is  shown  that  some  of  the  equations  usually  used  to  describe  Faraday 
rotation  are  Incorrect  for  the  case  of  propagation  through  an  tmlsotropic 
magnetic  field  when  the  cosine  of  the  angle  between  the  magnetic  field  and 
the  wave  normal  changes  sign  along  the  ray  path.  Such  a  case  Is  common  when 
transmissions  are  from  a  high  altitude  satellite  In  the  region  of  the  Ionosphere 
where  the  direction  of  propagation  is  newly  perpendicular  to  the  earth's 
magnetic  field.  Use  of  the  uncorrected  equations  can  result  In  misinterpretation 
of  the  data.  The  corrected  equations  are  given  and  are  shown  to  be  consistent 
with  measured  Faraday  rotation.  With  the  corrected  equations  used  In  a  ray 
tracing  program,  an  analysis  Is  made  of  errors  resulting  from  the  use  of  several 
approximations  near  the  transverse  region.  The  approximations  Investigated  are: 
(1)  quasl-longltxidlnal  propagation,  (2)  straight  line  propagation,  and  (3)  use 
of  an  effective  Ionospheric  height. 
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THTRODUCTIOH 

During  the  suamer  of  1962,  recordings  were  made  of  the  satellite  radio 
signals  from  the  Transit  ItA  satellite  at  a  low  latitude  receiving  station. 
During  the  analysis  of  the  Faraday  rotation  of  these  signals.  It  was  found 
that  sample  calcxilatlons  with  a  ray  tracing  program  gave  results  which  were 
not  always  consistent  with  the  measured  data.  It  was  found  that  whenever  the 
earth’s  magnetic  field  Is  perpendicular  to  the  ray  path  at  some  point  between 
the  transmitter  and  receiver,  the  equations  generally  iised  to  describe  the 
Faraday  rotation  are  Incorrect.  In  Section  I  of  this  report  the  correct 
equations  are  derived,  and  In  Section  II  they  are  applied  to  analyze  errors 
resulting  from  the  use  of  severeil  different  approximations. 

I.  B4UATI0NS  OF  FARADAY  ROTATION 

Whenever  a  receiving  station  is  at  a  sufficiently  low  geomagnetic 
latitude  (less  than  about  1*0°),  there  is  a  portion  of  the  sky  where  radio 
signals  from  a  satellite  propagate  perpendicular  to  the  earth's  magnetic 
field  at  some  point  cd.ong  the  ray  path.  This  portion  of  the  sky  Is  referred 
to  In  this  paper  as  the  "transverse  region",  and  the  point  of  perpendicularity 
along  the  ray  path  as  the  "transverse  point".  When  the  satellite  Is  on  one 
side  of  the  transverse  region,  there  Is  a  conqonent  of  the  earth's  magnetic 
field  parallel  to  the  direction  of  propeigation,  as  In  the  case  of  vertically 
downward  propagation  in  the  northern  magnetic  hemisphere.  When  the  satellite 
Is  on  the  other  side  of  the  transverse  region,  the  magnetic  field  has  a 
component  cmtlparallel  to  the  direction  of  propagation  (parallel,  but  In  the 
opposite  direction),  as  In  the  case  of  vertically  downward  propagation  In  the 
southern  magnetic  hemisphere. 

To  understand  propagation  near  the  transverse  region,  It  Is  necessary  to 
consider  the  polarizations  of  the  "ordinary"  and  "extrsujrdlnary"  magnetolonlc 
modes  aiid  the  resultant  wave.  Ry  definition,  the  ordinary  ray  Is  that  which 
Is  unaffected  by  the  magnetic  field  when  propagation  Is  transverse  to  the 
field.  When  there  Is  a  longitudinal  component  of  the  field,  collision  effects 
are  negligible,  and  the  wave  frequency  Is  higher  than  the  gyromagnetlc 
frequency,  the  ordinary  wave  Is  less  affected  by  the  magnetic  field  than  Is 
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the  extraordinary  vave.  For  purely  Ipn^^tudlnal  propagation,  the  tvo  waves 
are  circularly  polarized  with  opposite  sense  of  rotation.  Referring  to 
Figure  lA,  the  direction  and  magnitude  of  the  force  felt  by  a  charge,  q, 
moving  perpendicular  to  the  magnetic  field  Is  given  by  the  equation 
F  =  q  V  X  Letting  y  *  y  "i  and  B  -  B  ^ ,  then  F^^  »  q(vj  Bj^  -  v^^  Bj)  or 
F  =  qvBl,so  the  charge  would  circle  counterclockwise  as  shown.  An 
electron  would  circle  In  the  opposite  direction,  (clockwise}  as  shown  in 
Figure  IB  (where  the  magnetic  field  Is  directed  Into  the  paper). 


Consider  two  clrc\ilarly  polarized  waves  of  opposite  sense  of  rotation 
propagating  parallel  to  the  magnetic  field.  Again  referring  to  Figure  IB, 
it  can  be  visualized  that  the  electrons  have  more  chance  to  be  Influenced  by 
a  clockwise  wave  than  by  a  counterclockwise  wave,  and  that  the  Interaction  Is 
a  maximum  when  the  wave  frequency  Is  equal  to  the  gyromagnetlc  frequency. 

Hence  the  clockwise  wave  Is  more  affected  by  the  magnetic  field,  and  therefore 
is  the  extraordinary  wave. 


To  determine  the  direction  of  Faraday  rotation,  let  \is  consider  the 
ordinary  and  extraordinary  waves  at  two  planes  which  are  perpendicular  to 
the  wave  normal  and  a  distance  di  apart  along  the  wave  normal.  This 
situation  Is  depicted  In  Figure  1C.  At  plane  A  at  time  t  ,  let  the  E  vectors 
of  the  two  waves  (E  and  E  )  be  In  the  dlarectlon  of  1.  Both  the  ordinary  and 

0  3C 

extraordinary  waves  are  circularly  polarized,  so  that  as  time  progresses,  E 
rotates  counterclockwise  at  the  wave  frequency,  I!  rotates  clockwise,  and  the 
resultant  E  Is  plane  polarized  In  the  plane  of  i  and  At  plane  B  at  time 
t^,  the  phase  of  E^  is  lagging  (clockwise)  by  an  angle  0^  and  E  is  lagging 
(counterclockwise)  by  an  angle  (Viewing  E^  at  time  t^,  it  Is  seen  that 

it  describes  a  right-hand  helix  in  space,  and  E  describes  a  left-hand  helix.) 
However,  since  the  wave  length  of  the  extraordinary  wave  Is  greater  than  that 
of  the  ordinary  wave,  0^^  <  The  resultant  !e  is  In  a  direction  making  em 


clockwise  from  1.  As  time  progresses,  E^  rotates 


angle  0  = 

counterclockwise  and  E^^  rotates  clockwise  at  the  wave  frequency  so  that  E 
Is  plane  polarized  with  the  polarization  plane  lagging  (clockwise)  by  em 
angle  0  from  its  position  at  plane  A.  Therefore,  the  plane  of  E  Is  twisted 
clockwise  In  space,  so  we  may  say  that  the  total  Faraday  rotation  Is 
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FIGURE  I.  FOR  EXPLANATION  SES  TEXT. 
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clockwise  when  propagation  Is  parallel  to  the  oagnetlc  field.  In  the 
opposite  case,  when  propagation  Is  antiparallel  to  the  magnetic  field, 
the  total  rotation  of  a  downcomlng  ray  la  counterclockwise. 

It  may  he  noted  that  In  the  case  of  parallel  propagation.  If  the  electron 
density  or  longitudinal  coi^ponent  of  the  magnetic  field  Is  decreasing  along 
the  ray  path,  there  will  he  decreasing  amounts  of  total  Faraday  rotation, 
so  the  plane  of  the  E  vector  at  the  ground  will  rotate  counterclockwise  as 
viewed  from  the  transmitter. 

Whenever  the  propagation  la  transverse  to  the  magnetic  field,  the  above 
analysis  Is  not  valid  since  the  ordinary  and  extraordinary  waves  are  no 
longer  circularly  polarized.  For  purely  transverse  propagation  they  are 
linearly  polarized  In  space  quadrat\ire.  The  resultant  wave  Is  then 
circularly  polarized  and  Faraday  rotation  cannot  he  defined  simply  as  the 
rotation  of  the  plane  of  polarization.  However,  In  the  Ionosphere,  only  In 
the  specled.  case  of  a  receiving  station  directly  on  the  magnetic  equator  can 
there  he  transverse  propagation  along  the  entire  ray  path  from  a  transmitter. 

In  the  more  general  case,  consider  a  receiving  station  in  the  northern 
hemisphere  at  a  sufficiently  low  latitude  that  transverse  propagation  can 
exist,  and  consider  a  transmitter  fixed  In  the  trsuisverse  region  of  the  sky. 
There  will  he  only  one  point  sdong  the  ray  path  where  propagation  Is  purely 
transverse.  Above  this  point  the  propagation  will  have  a  component  peurallel 
to  the  magnetic  field,  while  below  the  point  there  is  a  component  anti¬ 
parallel  to  the  field,  (in  the  actual  case  described  later,  the 
receiving  station  was  at  16%  geomagnetic  latitude  and  a  ^  me  transmitter 
was  at  an  sd-tltude  of  970  km.  The  transverse  point  was  at  an  altitude  of 
360  km.  At  93^  km  the  angle  between  the  ray  direction  and  the  magnetic  field 
was  83.2°  at  360  km  the  angle  was  90°,  and  at  109  Ion  the  angle  wm  92.9°.) 

Flg\ire  ID  depicts  the  magnetic  field  geometry  In  such  a  case  where  a  linearly 
polarized  signal  Is  transmitted  from  T.  At  point  D  there  Is  a  parsdlel 
component  of  the  field,  at  point  E  the  field  Is  transverse,  and  at  point  F 
there  Is  an  antlpcuallel  component.  At  point  D  the  Faraday  rotation  Is  building 
up  clockwise  and  the  polarization  of  the  slgnsQ.  Is  newly  linear.  Approwhlng 
E,  the  polarization  of  the  signal  changes  from  elliptical  to  circular,  while 
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the  Faraday  rotation  (determined  hy  the  position  of  the  major  axis  of  the 
ellipse)  continues  to  build  up  clockwise.  E  is  but  a  point  in  space  and  not 
a  region,  so  no  Fasraday  rotation  is  built  up  there.  Just  below  B  there  must 
be  an  elliptical  conqponent  to  the  polarization.  The  phase  of  each  characteristic 
wave  is  continuously  varying  across  E,  so  there  can  be  no  discontinuity  in  the 
Faraday  rotation  at  E.  However,  there  is  a  reversal  in  the  direction  of  the 
longitudinal  con^nent  of  the  magnetic  field  beyond  E,  l.e.,  there  is  an 
antiparallel  component  of  the  field.  Hence  the  ordinary  wave  describes  a 
left»hand  helix  (actually  elliptical  rather  than  circular),  the  extraordinary 
wave  describes  a  right-hand  helix,  and  thi.  Faraday  rotation  btiilds  up 
counterclockwise.  Therefore,  the  total  rotation  of  the  wave  at  the  receiver 
is  the  sum  of  the  clockwise  rotation  from  T  to  E  and  the  counterclockwise 
rotation  from  E  to  F.  This  stun  can  be  a  net  clockwise  rotation,  zero,  or  a 
net  counterclockwise  rotation. 

The  Faraday  rotation  is  usually  described  by  the  equation: 

“  -  5  (fo  - 

where  P  and  P„  are  the  phase  path  lengths  (in  wave  lengths)  of  the  ordinary 
and  extraordliiary  waves  respectively  and  0  is  the  Faraday  rotation  in  cycles; 

^o  “  X  Js^  *^0  ^°o'  “  X  ^ 

where  X  is  the  free  space  wave  length  of  the  transmitted  radiation,  |jl  and 
)ji^  are  the  indices  of  refraction  of  the  ordinary  and  extraordinary  waves 
respectively,  and  the  integrals  are  along  the  ray  paths  from  the  transmitter 
to  the  receiver.  Whenever  the  transmitted  frequency  is  higher  than  the 
plasma  frequency  and  collision  effects  are  negligible,  the  plus  sign  in  the 
Appleton-Hartree  equation  is  used  for  the  ordinary  wave  and  the  minus  si£pi 
for  the  extraordinary  wave.  In  this  cMe  is  edways  greater  than  or  equal 
to  so  that  P^  >  Pjj  and  by  Equation  (l),  n  >  0.  Equation  (2)  does  not 
indicate  whether  the  waves  describe  right  or  left-hand  helices  and  give  only 
the  absolute  number  of  wave  lengths.  Hence  Equation  (l)  gives  the  sum  of  the 
absolute  values  of  the  clockwise  and  counterclockwise  rotation,  which  is  not 
indicative  of  the  total  rotation  between  the  transmitter  and  receiver. 
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The  deficiency  In  Equations  (1)  and  (2)  can  be  corrected  by  defining 
the  spatial  rotation  in  one  direction,  aay  counterclockwise,  to  be  positive. 

With  this  definition,  when  a  component  of  the  magnetic  field  is  parallel  to 
the  ray  path  (as  in  vertically  downward  propagation  in  the  northern  hemisphere), 
the  ordinary  wave  is  negative  (clockwise),  the  extraordinary  wave  is  positive 
(counterclockwise),  and  the  total  rotation  is  negative  (clockwise).  When 
there  is  an  antiparallel  component  of  the  magnetic  field,  all  three  are 
reversed. 

Also  Equations  (l)  and  (2)  then  become: 


where  8  =  +  1  when  the  ray  path  is  parallel  to  the  longitudinal  component ' of 
the  magnetic  field  and  8  =  -  1  when  the  ray  path  is  antiparallel  to  the  field. 

In  the  case  illustrated  by  Figure  ID,  8  =  +  1  from  T  to  B  and  8  =  -  1  from 
B  to  R. 

In  practice,  it  would  be  difficult  to  do  ray  tracing  with  the  use  of 
Equations  (3)  and  (4).  This  is  because  it  would  be  necessary  to  determine 
the  transverse  point  within  one  wave  length  of  the  transmit ied  radiation  on 
both  the  ordinary  and  extraordinary  ray  paths  in  order  for  the  resulting 
Faraday  rotation  to  be  accurate  to  within  one  cycle.  In  a  ray  tracing 
program  such  as  described  by  Lawrence  and  Posakony  (19^1)  or  Dulk  and  Dean  (1962), 
this  would  require  layer  thicknesses  nesur  the  transverse  point  of  about  6  meters 
for  a  me.  signal.  Round-off  errors  for  these  small  layers,  eis  well  as  the 
possible  non-convergence  of  the  iteration  process  for  computing  refraction, 
would  make  the  procedure  Impractical.  Therefore,  another  scheme  has  been 
devised  for  ray  tracing  in  the  transverse  region.  Advantage  is  taken  of  the 
facts  that  magnetoionlc  splitting  is  a  minimum  for  nearly  transverse 
propagation  and  that  the  "no  magnetic  field"  ray  path  is  between  the  ordinary 
ray  path  and  the  extraordinary  ray  path,  so  that  least  error  is  Introduced 
by  tracing  along  the  "no  magnetic  field"  ray  path.  The  Faraday  rotation  is 
computed  with  the  equation: 
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(5) 


"■fe 

where  refraction  but  not  magnetolonlc  splitting  Is  taken  Into  account  along 
the  "no  magnetic  field"  ray  path.  This  equation  may  be  used  whenever  the 
transverse  point  occurs  somewhere  along  the  ray  path.  At  the  transverse 
point,  6  changes  from  +1  to  -1.  Using  Equation  (5),  the  transverse  point 
need  not  be  determined  as  closely  as  when  using  Equations  (?)  and  (k),  and 
the  sign  of  U  Indicates  whether  the  total  rotation  Is  counterclockwise  or 
clockwise. 

There  are  two  special  cases  of  Interest.  For  the  first,  the  treuisverse 
point  coincides  with  the  bottom  of  the  Ionosphere.  In  this  case  the  limiting 
polarization  of  the  wave  Is  circular,  and  this  should  be  noticeable  on  signal 
strength  records  as  a  decrease  In  depth  of  the  Faraday  nulls.  This  occurrence 
usually  will  not  be  coincident  with  the  true  zero  total  Faraday,  which  occurs 
when  the  transverse  point  occurs  near  the  middle  of  the  ionosphere.  For  the 
second  case,  the  transverse  point  Is  at  the  transmitter  height  and  the  plane 
of  polarization  of  the  transmitted  wave  is  In  the  plane  of  the  magnetic 
meridian.  Then  only  the  ordinary  mode  will  be  transmitted.  At  the  receiver 
the  signal  will  be  almost  clrcuOarly  polarized.  Should  the  plane  of  polarization 
of  the  transmitted  wave  be  perpendicular  to  the  magnetic  meridian  at  the 
transmitter,  only  the  extraordinary  mode  would  be  transmitted  and  a  circularly 
polarized  signal  of  opposite  sense  of  rotation  would  be  received. 

H.  ANALYSIS  OF  ERROR  INVOLVED  IN  VARIOUS  APIROXIMATIONS 

A  thorotigli  Investigation  was  made  of  one  satellite  pass  during  which 
transverse  propagation  occurred.  The  primary  purpose  of  the  Investigation 
was  to  determine  the  magnitude  of  the  errors  Involved  In  making  various  - 
approximations.  The  pass  investigated  was  revolution  5183  of  the  Transit 
I4A  satellite  which  was  transmitting  on  5^  me  and  324  me.  The  satellite 
signals  were  received  at  a  low  latitude  receiving  station.  (Johnston  Island, 

16°  N,  190°  E;  geomagnetic  latitude  15°  N. )  This  pass  was  from  north  to  south, 
the  closest  approach  of  the  satellite  being  about  500  km  east  of  the  receiver 
at  10:46  AM  Johnston  Island  standard  time  and  the  satellite  altitude  was  about 
960  km.  During  the  pass  the  Ionosphere  was  relatively  quiet,  the  geomagnetic 
planetary  A  Index  was  10. 
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FIGURE  2 

SUBSATCLLITE  ELECTRON  CONTENT 
AND  EFFECTIVE  IONOSPHERIC  HEIGHT 
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The  ray  tracing  progzam  used  In  the  analysis  vas  that  of  Dulk  and 

Dean  (1962).  Thou^  an  option  for  ray  tracing  through  an  Ionosphere  varying 

In  three  dimensions  Is  available  In  the  program,  coiqiutatlonB  are  slogtllfled 

when  the  Ionosphere  model  Is  spherically  stratified.  Revolution  ^183 

exhibited  a  nearly  constant  subsatelUte  electron  content  over  a  considerable 

portion  of  the  pass,  enabling  vise  of  a  spherlcedly  stratified  lonospdiBre 

model  In  the  ray  tracing  program.  A  preliminary  Investigation  using  the 

hybrid  dispersive  Doppler-Faraday  rotation  method  of  De  Mendonca  (196S) 

17 

Indicated  that  the  sxibsatelllte  electron  content  vas  about  2.2  x  10 
electrons/meter  and  that  the  electron  content  Increased  sli^tly  to  the 
soxith  of  the  receiving  station.  Figure  2  shows  the  electron  content  poroflle 
as  fo\md  by  the  hybrid  method. 

The  electron  density  profile  used  vas  In  two  sections: 

1.  Up  to  (the  height  of  f2^^)  the  profile  vas  derived  from 

an  lonogram  taken  by  a  low  latitude  lonosonde  (Puerto  Rico)  under  similar 
Ionospheric  conditions. 

2.  Above  the  profile  vas  eissumed  to  be  that  of  a  Chapman 

layer  with  a  scale  height  at  of  42  km  (indicated  by  the  layer  half 
thickness  as  shown  on  the  lonogram)  and  a  scede  height  gradient  of  .12  km/km. 
The  Chapman  scale  height  and  sceLLe  height  gradient  were  discussed  In  detail 
by  Marks,  et  al,  (1963),  and  were  shown  to  be  realistic. 

For  the  purposes  of  the  major  pert  of  the  present  Investigation,  the 

model  Ionosphere  profile  was  not  required  to  be  In  exact  agreement  with  the 
profile  actually  existing  at  the  time  the  satellite  was  tracked.  It  was 
necessary  only  that  the  same  model  be  used  In  testing  e6u:h  of  the  various 
approximations.  However,  to  enable  cos^Artson  with  meastved  data.  It  was 
desirable  for  the  model  profile  and  the  actual  profile  to  be  similar  In  shape, 
and  to  contain  the  same  subsatellite  electron  content.  Therefore,  the  electron 
density  model  was  nozmallzed  to  give  the  subsatellite  electron  content  of 
2.2  X  10  '  electron/meter  to  be  In  agreement  with  the  hybrid  measurement. 

Also  the  shape  of  the  model  profile  and  electron  density  at  were  later 

found  to  be  In  agreement  with  lonosonde  data  taken  at  Johnston  Island  during 

the  satellite  pews. 


FIGURE  3 
FARADAY  ROTATION 
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The  first  step  In  the  Investigation  involved  ray  tracing  the  ^  me 
signal  from  satellite  to  receiver  at  about  20  satellite  locations •  These 
locations  were  distributed  from  horizon  to  horlzon>  but  the  majority  were 
in  the  transverse  portion  of  the  sky.  In  Figure  3f  the  Faraday  rotation  as 
consiuted  In  the  ray  tracing  Is  compared  with  the  measured  Faraday  rotation. 
Because  the  time  of  zero  meeisured  Faraday  rotation  cannot  usually  be  determined 
from  the  record.  It  was  arbitrarily  chosen  to  be  at  206  seconds,  the  time  of 
zero  cosq^uted  Faraday.  It  la  seen  in  the  Figure  that  there  Is  good  agreement 
between  the  two  sets  of  data  over  most  of  the  IntervcG.  of  measurements.  Only 
at  later  times,  with  the  satellite  to  the  south,  does  the  measured  Faraday 
become  greater  than  the  Faraday  cosiputed  from  the  model,  indicating  increasing 
subsatellite  electron  contents  to  the  south. 

The  magnetic  field  model  used  in  the  comgputatlons  was  the  6th  order 
spherical  harmonic  description  using  coefficients  of  Jensen  and  Cain  (1962) 
for  Epoch  i960.  Using  this  model.  It  was  found  that  when  signals  were  first 
received  at  the  beginning  of  the  pass  (089  sec.),  there  was  a  transverse  point 
along  the  ray  path  In  the  Ionosphere  at  an  altitude  of  9^  As  the 
satellite  traveled  southward  In  its  orbit,  a  transverse  point  continued  to 
exist  somewhere  along  the  ray  path  for  a  period  of  ll(0  sec,  during  which  the 
satellite  traveled  from  29.2°  N,  178.1v°  B  to  21.7°  N,  181.7°  E. 

Figure  4  iUustrates  the  error  involved  in  Ignoring  refraction  of  the 
ray.  0  was  obtained  by  using  Equations  (3),  (4),  and  (^),  the  same 

equations  used  to  obtain  n  Trace’  however,  the  Integral  was  taken  along 
the  stredght  line  path  from  transmitter  to  receiver  Instead  of  along  the 
refracted  ray  paths.  In  all  but  the  transverse  region,  the  resulting  error 
Is  seen  in  Figure  4  to  be  less  than  about  2^,  with  the  error  increasing  at 
later  times  with  the  satellite  approaching  the  horizon.  However,  near  the 
point  where  0  =  0,  which  occurs  with  the  transverse  point  at  an  altitude  of 
about  3^  4m,  Ignoring  refraction  causes  the  error  to  be  greater  than 
for  a  period  of  about  30  seconds.  During  the  early  portion  of  the  perns,  with 
the  transverse  point  between  ItOO  km  and  96O  km,  the  error  Is  about  6^.  The 
reason  for  the  sensitivity  to  refraction  effects  In  the  transverse  region  Is 
that  the  sllf^t  difference  between  the  true  ray  direction  and  the  straight 
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ilna  path  causes  the  treuusverse  point  to  be  shifted  upward  or  downward  to  a 
considerably  greater  extent.  The  computed  Faraday  rotation  Is  therefore 
altered  by  a  disproportionate  amount. 

The  error  Involved  in  \ising  the  quasl-longitudloal  appraodoatlon 
throughout  the  transverse  region  is  Illustrated  in  Figure  3*  The  equation 
utilizing  the  quasi-longitudinal  approximation  which  was  used  to  conqnite 
the  Faraday  rotation  Is 

where  K  is  a  constant  ^  f  is  the  transmitted  frequency,  H  is  the  magnetic  field 
strength,  d  is  the  angle  between  the  sagnetie  field  and  the  wave  normal,  N  Is 
the  electron  density,  and  the  integral  is  eiLong  the  "no  magnetic  field"  ray  path 
from  transmitter  to  receiver.  The  error  of  1^,  which  is  seen  In  Figure  3  to  be 
nearly  constant  over  much  of  the  pass,  can  be  attributed  to  use  of  a  constant  K 
which  was  1^  too  high.  Only  in  the  ten  second  interveil  about  0  -  0  is  the  error 
greater  than  10^.  Kven  in  the  early  portion  of  the  satellite  pass,  with  the 
transverse  point  occvxrring  somewhere  along  the  ray  path,  use  of  the  Q.L. 
approximation  resvilts  in  little  error.  Conqierlng  Figure  4  with  Figure  It 
Is  seen  that  more  error  is  incurred  in  Ignoring  refraction  than  in  using  the 
Q.L.  apprcKlmatlon,  and  the  interveil  of  large  error  is  longer. 

Figure  6  Illustrates  the  error  in  vising  the  Q.L.  approximation  along  the 
strcd.ght  line  path  from  transmitter  to  receiver.  The  result  is  essentially  a 
superposition  of  the  errors  in  Figure  4  euad  Figure 

Wien  using  the  Q.L.  approximation  above,  the  factor  H  cos  9  was  retained 
within  the  integral  sign.  Often  It  is  desirable  to  find  a  mean  value  of  the 
factor  H  cos  9  (or  H  cos  9  sec  X)  and  remove  it  from  under  the  integral  sign. 
Then  Equation  (6)  Is  changed  to 


E  cos  0  N  ds 


(6) 


n 


(H  cos  9  sec  X)_ 

f2 


(7) 


Where  is  the  proper  mean  height  of  the  ionosp^re,  X  is  the  angle  between 

the  ray  path  and  the  zenith  at  height  H. ,  and  f  dh  is  the  subsatellite 

1  J  o 
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electron  content.  Figure  7  Ulustratea  the  error  resulting  from  use  of 
Equation  7  when  the  mean  or  effective  Ionospheric  height  Is  chosen  to  he 

330  km.  Again,  over  much  of  the  pass  the  error  Is  less  than  2$,  Errors  of 
over  10^  occur  during  an  Interval  of  about  20  seconds  near  the  time  uhen 
n  0.  During  early  tines  the  error  Is  seen  to  be  about  6$,  with  Equation  7 
overestimating  the  rotation. 

The  resvilts  of  assuming  the  mean  Ionospheric  hel^t  to  be  too  low,  at 
300  km,  are  shown  In  Figure  8.  Here  the  error  Is  seen  to  be  about  over 
much  of  the  peiss,  with  Equation  7  ovearestlmatlng  the  amount  of  Faraday 
rotation.  Near  the  tine  when  n  =  0,  there  Is  an  Interval  of  about  73  seconds 
when  the  error  Is  greater  than  UO^.  During  early  tines.  Equation  7 
\mdereatlmated  the  amount  of  rotation  by  about  10^.  Comparison  of  Figure  8 
with  Figure  7  shows  that  a  mean  Ionospheric  height  of  about  3^  km  would  be 
the  best  for  the  particular  Ionospheric  structure  used  In  the  computations. 

The  results  of  overestimating  the  mean  Ionospheric  height  by  setting 
=  400  km  and  using  Equation  7  are  shown  In  Flginre  9«  Again  there  Is  an 
error  of  about  over  most  of  the  pass,  with  the  rotation  ceQ.culated  with 
Equation  7  being  too  low.  Near  the  time  when  0=0,  there  Is  an  Interval  of 
about  25  seconds  when  the  error  Is  greater  then  From  the  beginning  of 

the  pass  until  213  seconds,  an  Interval  of  about  123  seconds,  the  error  Is 
greater  than  20^.  At  early  times.  It  Is  seen  that  Equation  7  overestimates 
the  amount  of  rotation. 

Figure  10  Is  a  composite  of  the  percent  errors  shown  In  the  preceding 
three  figures.  It  Is  notable  that  there  Is  a  point  at  about  263  seconds 
where  the  Faraday  rotation  given  by  Equation  7  Is  Independent  of  the  altitude 
chosen  for  H^.  This  la  the  time  when  the  feu^tor  H  cos  9  sec  X  changes  from 
a  decreasing  function  of  altitude  to  an  Increasing  ftinctlon  of  altitude. 

Such  a  point  has  been  found  for  all  satellite  pcwses  so  far  Investigated. 
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CONCmSIOH 


It  has  1)6611  shown  that  the  usual  equations  for  decrlblng  Faraday 
rotation  are  Incorrect  whenever  there  Is  a  change  of  sign  of  cos  9  somewhere 
edong  the  ray  path*  The  corrected  equations  are  given,  and  are  shown  to  be 
consistent  with  measured  data.  Use  of  the  uncorrected  equations  can  lead 
to  In^oper  Interpretation  of  Faraday  rotation  near  the  transverse  region. 

With  the  corrected  equations  used  In  a  re^  tracing  program,  em  analysis 
of  errors  Involved  In  several  approximations  shows  that  Ignoring  refraction 
can  result  In  considerable  error  near  the  trtuisverse  region  even  at  a  frequenay 
as  high  as  me,  that  use  of  the  quasi-longitudinal  apptroxlmatlon 
results  In  less  error  then  Ignoring  refraction,  and  that  the  choice  of  the 
proper  mean  Ionospheric  height  Is  quite  Important  when  deducing  Ionospheric 
electron  content  from  measurements  of  Faraday  rotation. 
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